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ABSTRACT 

Minor accretion events with mass ratio -&f 8 at : -^hoBt — 1:10 are common in the 
context of ACDM cosmology. We use high-resolution simulations of Galaxy-analogue 
systems to show that these mergers can dynamically eject disc stars into a diffuse 
light component that resembles a stellar halo both spatially and kinematically. For 
a variety of initial orbital configurations, we find that ~ 3 — 5 x 10 s Mq of primary 
stellar disc material is ejected to a distance larger than 5 kpc above the galactic plane. 
This ejected contribution is similar to the mass contributed by the tidal disruption of 
the satellite galaxy itself, though it is less extended. If we restrict our analysis to the 
approximate solar neighborhood in the disc plane, we find that ~ 1% of the initial disc 
stars in that region would be classified kinematically as halo stars. Our results suggest 
that the inner parts of galactic stellar halos contain ancient disc stars and that these 
stars may have been liberated in the very same events that delivered material to the 
outer stellar halo. 
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1 INTRODUCTION 

The canonical paradigm describing the emergence of 
outer stellar halos around galaxies involves the accretion 
and tidal disruption of satellite systems llSearle &; Zinn 



19781; iJohnston et al.lll996l; ICote et"ai] l200d:lBullock et al 
2001al: iBullock fc Johnston! 120051: l Abadi et all l200fj; 



De Lucia fc Helmill2008l ; I Johnston et alj|2008l ; ICooper et ail 
20091 ). 



Empirical evidence suggests that these pre- 
dicted mergers have indeed populated the outskirts of 
the Milky Way and M31 halos, with diffuse, stream- 
like signatures apparent at low surface bright ness (r e- 

l2009al lbl: 
references 



cently, iMcConnachie et al.l 12009: 


Gilbert et al.l 


1 St arkenburg et al.l 


2009: Bell et al. 


20081. and 



therein). 

Perhaps unsurprisingly, significant diffuse stellar halo 
components also appear to be pervasive in the population of 
nearby disc galaxies outsid e the Local Group (ISackett et alJ 
19941 : Morrison et al.lll997l; [Lequeux et ai1ll998l : lAbe et all 
19991 : Izibetti fc Ferguson! l2004h. The GHOSTS survey de- 



scribed by de Jong et al. ( 2007a! ) has established the pres- 
ence of relatively luminous and extended halos surrounding 
the massive disc galaxies NGC253, NGC891, and M94; the 
1^-band surface brightness fiv of this material is typically 
~ 28 — 29 mag arcsec -2 at large heights (~ 20 — 30 kpc) 
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above the disc plane. This quantity of diffuse light can- 
not easily be reconciled with a model where most of the 
inner, richer stellar halo is built by the destruction of dwarf 
satellite galaxies accreted at early times, since these old 
and faint subhalos can add at most only a few percent 
to the total luminosity of the host galaxy and produce 
metal-poor stellar halos jPurcell et al.ll200l[20o3 ). This ac- 
creted light is typically deposited over such a large vol- 
ume that the resulting ancient halo's surface brightness 
is far below that observed in the less-extended spheroids 
enveloping many massi ve disc galaxies. Faint st ellar halo 
material around M33 (IMcConnachie et al.l I200R see also 
IChapman et al.ll2006l : iKalirai et al.l l2006l : iHelmil |2008h and 
diffuse extended ligh t around dwarf galaxies like the LMC 
ijMinniti et all 120031 ) ar e even more difficu lt to explain in 
the accretion scenario jPurcell et alj|2007h . These general 
argument s, along with direct eviden ce that the inner halo 
of M31 dGuhathakurta et all 120051 ) and the Milky Way 
jCarollo et . al.ll2007h are populated by at least two distinct 
components, make it clear that a multivalent picture of inner 
stellar halo formation must emerge. 

There are many sources for hot stellar halo mate- 
rial besides direct accretion. Gas-rich mergers in the early 
epochs of galaxy assembly are likely important for the cre- 
ation of inner halo material (|Brook et all l2004al lbh. An- 
other related possibility is that in-situ stars which form at 
high redshift (z > 3) become ejected to large galactocen- 
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Figure 1. Surface brightness maps of our primary system, viewed edge-on: the initial model is visualized in the left panel, while the 
center panel shows only stars belonging to the primary galaxy following the prograde infall with orbital inclination of 60° , and the right 
panel subsamples the accreted stars only. In all cases, surface densities have been converted into surface brightnesses using a stellar- 
mass-to-light ratio M*/L = 3; note that the brightness limits have been chosen in order to accentuate stellar halo substructure at the 
expense of saturating the galactic disc region. 



trie radii by subseq uent major mergers between 2 < z < 3 
l|Zolotov et all 12009). In small galaxies, supernovae can af- 
fect the galactic structure in such a way that primary 
galaxy stars become liberated int o an extended, stellar halo 
distribution (IStinson et alj 120091 ) . Here, we focus on mi- 
nor mergers of the type that shoul d be common in th e 
last ~ 10 Gyr of galactic assembly l|Stewart et al.| [2008). 
These accretion events involve satellites massive enough to 
severely heat the host's stellar distribution, but not massive 
enough to compl etely eliminate the presence of the primary 
disc component jKazantzidis et alj 20081; iRead et af1l2008l; 
Villalobos fc Helmill2008l; IPurcell et aLll2009al : lMoster etaU 
20091 ; [kazantzidis ct al. 20j)|). 



Using a suite of high-resolution collisionless simulations 
as well as resimulated test cases involving a full treatment 
of hydrodynamical physics, we investigate the process of dy- 
namical stellar ejection induced by the accretion of satellites 
onto a thin galactic disc, and demonstrate that disc stars can 
become members of the stellar halo, as it is defined both spa- 
tially and kinematically, via local and global heating modes. 
In order to determine the dynamical response of a disc to 
a cosmologically-motivated accretion event, we simulate the 
interaction between a Galaxy-analogue primary system with 
virial mass Mhost — 1O 12 M0 and an infalling satellite sub- 
halo much more massive than the Galactic disc itself, M sa t ~ 
10 U M S ~ 3M disc . This mass ratio M sa t:M hoat ~ 1:10 has 
been sh own, via analy tic formulations jPurcell et alj [20071 . 
see also Zentncr 200 7|) as well as cosmological s imulations 
jStewart et al.ll200Sl ; iBovlan-Kolchin et al.ll2009l ). to repre- 
sent the dominant mode of mass delivery into CDM halos, 
and is therefore of immediate co ncern for the questio n of 
disc survivability (as addressed bv lPurcell et al.l l2009a). 

The majority of the simula tions we pre s ent he re 
have been discussed elsew here (IPurcell et al.l l2009al lbl: 
iKazantzidis et alj 120091 ) . In IPurcell et alj (|2009al ) we ar- 
gued that the resultant discs are likely not good analogues 
to the Milky Way. Nevertheless, it very well may be that 
the Milky Way is unusually thin, cold, and old. The re- 
sultant discs are not necessarily inconsistent with the ma- 
jority of external disc galaxies (|Moster et al.l [20091 ) . which 
do seem to have thicker scale heights than the Milky Way 



jYoachim fc Dalcantonll2006l ). though observational uncer- 
tainties associated with dust-lane obscuration and other ef- 
fects are important caveats for these considerations. Our 
goal here is to examine a plausible scenario for the emer- 
gence of an inner stellar halo component in most galaxies. 
They should not be compared directly to the Milky Way 
without some caution (see §3.1). In §2, we outline our nu- 
merical experiments, following with a presentation of our 
results in §3, reserving §4 for conclusions and discussion. 



2 METHODS 

We utilize high-resolution, multi-million-particle simulations 
in order to analyze the response of a thin stellar disc to the 
infall and accretion of a satellite galaxy one-tenth as massive 
as the host halo itself. Using fully self-consistent numerical 
models for both the primary and satellite galaxies, as well 
as initial conditions for accreting subhalos drawn from dis- 
tributions defined by large-scale cosmological simulations, 
we investigate a range of orbital infall inclinations. Alterna- 
tive descriptions of our simulations are provided in Purcell 
et al. (2009a, b); we now review the essential features for 
completeness. 

All of our numerical investigations use the multi- 
stepp ing, parallel, tree JV-body code PKDGRAV (|Stadell 
2001), in which we set the gravitational softening length to 
e = 100 pc and 50 pc for dark matter and stellar particles, 
respectively. Each simulation is collisionless, with the ex- 
ception of two experiments in which we ascertain the effect 
of a full treatment of hydrodynamical physics on the ejec- 
tion of both old and newly-formed disc stars from a primary 
galaxy with a modest initial gas fraction f g = 0.1. The latter 
simulation w as performed with t he TreeSPH iV-body code 
GASOLINE <|Wadslev et al.ll2004) and is part of a campaign 
designed to investigate the effect of a dissipative component 
on the dynamical response of thin galactic discs subject to 
bombardment by halo substructure (jKazantzidis et al ]|2010l . 
in preparation). In these cases, the gravitational softening 
length for the gas particles was set to e = 50 pc. 

Our primary galaxy components, drawn from the 
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Figure 2. Central surface brightness profiles as measured along 
the minor axis for each simulation endstate, in both the colli- 
sionless (thick lines) and hydrodynamical (thin lines) regimes, as 
a function of orbital inclination angle. Here, as in Figure [T] we 
assign a surface brightness /i based on the stellar-mass-to-light 
conversion factor M* / L = 3. For reference, at heights z > 5 kpc, 
we plot a profile with arbitrary normalization and power-law in- 
dex n = —2.5. 



fully self-consistent distribution functions advanced by 
IWidrow et alj |2008l ). represent equilibrium solutions to the 
coupled collisionlcss Boltzmann and Poisson equations, and 
are thus ideally suited for simulating the complex dynam- 
ics involved in the accretion of a massive satellite galaxy. 
The host system consists of a massive disc as well as a 
central bulge following a Sersic profile of effective radius 
R e = 0.58 kpc and index n = 1.118 (with mass Afbuigc = 
9.5 x 10 9 Mq distribu ted among 5 x 10 5 pa rticles); this is the 
model labeled Gl in lPurcell et all (1200981 ). The stellar disc 
is initialized with an exponential scale length Rd = 2.84 kpc 
and a vertical distribution described by a sech 2 function 
with scale height Zd = 0.43 kpc, and contains 10 6 particles 
comprising a total disc mass Mdi sc = 3.6 x 1O 1O M0. 

In the hydrodynamical runs, a fraction f g = 0.15 of 
the stellar disc mass is converted into gas, and this gaseous 
component is constructed with the same initial density dis- 
tribution as the stellar disc; the adopted method ology will 
be described in detail in iKazantzidis et al.l (|2010l . in prepa- 
ration). We include atomic cooling for a primordial mixture 
of hydrogen and he lium, a nd th e star formation algorithm 
is based on that of iKatj (| 19921 ). in which gas particles in 
cold and dense Jeans-unstable regions as well as convergent 
flows spawn star particles at a rate proportional to the local 
dynamical time. In our particular application, gas particles 
are eligible to form stars if their density exceeds 0.1 cm -3 
and their temperature drops below 1.5 x 10 4 K. Feedback 
from supe rnovae is treat e d usin g the blast-wave model de- 
scribed in Sti nson et al.l l|2006l ); in this model, the energy 
deposited by a Type-II supernova into the surrounding gas 
is 4 x 10 50 erg. This choice of parameters and numerical 
techniq ues produces realistic g alaxies in cosmological simu- 
lations (|Governato et al1l2007 ). Each primary galaxy model 



is embedded in a dark host halo composed of 4 x 10 par- 
ticles and arranged in a structure fo l lowing the canonical 
NFW density profile of lNavarro et al.l |l996), with scale ra- 
dius r s = 14.4 kpc and virial mass Mhost — 1O 12 M0. We 
note that our primary galaxy model is a Ga lactic-analogue 
system chosen from the lWidrow et akl l|2008l ) set of models 
with characteristics corresponding closely to criteria deter- 
mined observationally for the Milky Way, and that our pre- 
cise choice of the initial parameter set minimizes the effects 
of secular evolution such as the formation of a strong central 
bar, as well as artificial heating induced by the interaction 
of disc particles with more massive halo particles. 

The satellite galaxy simulated in each isolated accre- 
ti on event is initialized using the self-consistent formalism 
of iKuiiken fc Dubinskil Jl995), with which we embed a to- 
tal stellar mass M* = 2.2 x lO 9 A/0 (adopting the value for 
M*/M sat at z ~ 0.5 derived in the nu mber-density matching 
exercise o f lConrov fc Wechslerllioogh . comprised of 10 s par- 
ticles populating a spheroidal distri bution with Sersic in dex 
n ~ 0.5 (according to the findings of lvan Zee et al.ll2004l . for 
dSph shape parameters as a function of magnitude in Virgo 
cluster member systems), into a dark subhalo composed of 
10 6 particles with a virial mass 1O 11 M0 and having a den- 
sity structure well-fit by an NFW profile with concentration 
c v i r ~ 14 at z ~ 0.5. This value of c v ir is consistent with 
subhalo distributions drawn f rom cosmological s imulations 
ijBullock et al.ll2001bl . see also lMaccio et alJl2007h . 

To investigate the dependence of heating efficiency dur- 
ing a satellite-disc interaction on the orbital parameters of 
the accreting subhalo, we simulate a range of inclination an- 
gles for each primary galaxy model: 9 — 0°, 30°, 60°, and 90° 
(where 9 is defined as the angle between the angular momen- 
tum axes of the galactic disc and the infall orbit). In each 
experiment, the satellite galaxy is initially placed at a large 
distance r ~ 120 kpc from the center of the host halo, to 
minimize the shock imparted to the primary galaxy parti- 
cles by the sudden appearance of this new potential well. 
The initial velocity vector for the subhalo is motivated by 
cosmological analysis of substructure accretion, in which the 
distributions of radial and tangential velocity components 
for satellite systems (v r and v t ) peak respe ctively at 90% 
and 60% of the host halo 's virial velocity l)Bensonl 120051 : 
iKhochfar fc Burker j 120061 ) ; therefore, our subhalo vectors 
have v r = 116 km/s and vt = 77 km/s. 

All fiducial simulations involve orbits that are prograde 
with respect to the primary galaxy's rotation; we also per- 
form a retrograde simulation for the inclination angle 9 — 
60° in order to assess the degree to which this case decreases 
the efficiency of disc heating while increasing the amount by 
which the primary galaxy is tilted away from its initial angle 
of repose, as expected according t o previous experiments in- 
volving satellite-d i sc int eractions l|Velazquez fc Whiteill999l ; 
IKazantzidis et all 120091 ) . Each prograde simulation is al- 
lowed to evolve for a total of 5 Gyr, by which time the 
subhalo has fully coalesced with the host halo and the heat- 
ing process is virtually complete, leaving the galaxy in a 
quasi-steady state; the elapsed timescale for the retrograde 
model is slightly longer at 7 Gyr, since dynamical friction 
is less efficient in this case. Lastly, the hydrodynamical sim- 
ulations employ prograde orbits with inclination angles for 
the infalling satellite of 9 = 0° and 60°. 
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Figure 3. Minor-axis surface brightness fi (upper panels) and cumulative stellar mass M(> z) (lower panels) contributed by heated disc 
stars (solid lines) and accreted satellite stars (dashed lines) at height z from the galactic plane. 



3 RESULTS 

Regardless of orbital inclination, each fiducial accretion 
event ejects a significant portion of stellar mass from the 
disc, as shown qualitatively in the visualizations of Figure Q] 
This material reaches heights above the galactic plane con- 
sistent with the extended faint spheroids observed in many 
nearby massive disc galaxies, most notably in the diffuse 
inner halo of M31 (|Guhathakurta et all 120051 ; iKoch et all 
2008). The disc stars still located at the mid-plane of our 
simulated disc remnants are much dynamically hotter than 
the initially cold and thin stellar system, yielding a distri- 
bution of velocities consistent with a two-component disc as 
well as a non-trivial fraction of stars that have kinematics 
similar to those found in the stellar halo of the Milky Way. 

Figure [2] presents the minor-axis surface-brightness pro- 
files of the initial disc galaxy (solid, black) and the resultant 
discs after each encounter, as measured along the minor 
axis extending from the galactic centefl We see that the 
post-merger disc surface-brightness profile generally begins 
to show a distinct transition to a diffuse power-law compo- 
nent at fj, ~ 28 mag/ arcsec 2 (for a stellar- mass-to- light ratio 
M+/L = 3). Shown for comparison is a power-law slope of 
/i tx z~ 2 ' 5 , which is indicative of that expected for a stellar 
halo. As measured and reported by IPurcell et al.l |2009a) 



We choose the central region since it is the brightest and thus 
most well-constrained area in observation of external galaxies. 
The minor-axis profile does not change in slope significantly as 
galactocentric radius increases, although the disc outskirts at 
-R> 10 kpc show pronounced flaring at faint surface brightnesses. 



using a two-component sech 2 fitted profile, the scale height 
of the remaining thin-disc component varies along the range 
2d ~ 1 — 2 kpc according to the satellite's orbital inclination. 
We note that the thicker component of the fit invariably has 
a very large scale height (^diffuse ~ 4 — 7 kpc), reflecting the 
large amount of energy imparted through global and direct 
heating modes during the accretion event. 

Generally, our experiments result in a significant por- 
tion of the initial disc's stellar mass being ejected by the 
satellite-disc interaction into a diffuse spheroidal compo- 
nent at heights greater than z > 5zd ~ 7 kpc. Figure 
shows the surface brightness profiles (as in Figure [2]) as 
well as the cumulative stellar mass found above a given 
height z from the disc plane, with both quantities having 
been split into the components contributed by the heated 
disc stars and the accreted satellite stars. We see that low- 
latitude accretion events (6 = 0°,30°) heat the discs so 
efficiently that a stellar mass equivalent to roughly 1% of 
the initial disc mass, 

M haio ^ 2 - 3 x 10 8 M Q , is ejected 
to heights larger than 7 kpc. For high-latitude accretion 
events (0 = 60°, 90°), tidally-stripped stars from the ac- 
creted satellite dominate the population, and disc-star ejec- 
tion contributes M^ al ° ~ 10 s M . Closer to the galactic mid- 
plane, the very thick stellar disc becomes the prevalent fac- 
tor regardless of inclination angle; approximately 1O 9 M0 in 
initial-disc stars alone are ejected to a height larger than 
z ~ 4 kpc, across the full radial extent of the disc. 

Despite some thin-disc regrowth via star formation in 
the hydrodynamical simulations, the vertical disc thicken- 
ing is even more pronounced in those cases than in the 
collisionless experiments, possibly because the initial stellar 
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Figure 4. The solar-neighborhood Toomre energy diagram for the initial model (left panel) and the endstate system including heated 
disc stars (center panel) and accreted stars (right panel) following the satellite galaxy's prograde infall at an orbital inclination of 0° 
onto a primary galaxy with an initial gas fraction f g = 0.15; the quantity T is the quadratic sum of the radial and vertical space velocity 
components U = vr and W = v z , and V is the rotational velocity v^. Disc stars are subsampled from the galactic discs in a very thin 
annulus (100 pc in width and 300 pc in height) centered vertically on the plane of the disc and sweeping through the circumference 
defined by the solar radius Rq = 8 kpc; the annulus width and height are each tripled for sampling of the accreted stars, to minimize 
numerical noise. The origin of the x-axis is set in each case by the local standard of rest, defined as the mean rotational velocity V of 
the stellar subsample. We define stars inside the enclosing dashed line with a radius of 70 km/s to be kinematic members of the thin 
disc, while stars outside this region and enclosed by the dotted line with a radius of 180 km/s are members of the thick disc (as in 
Nissen fc Schuster! l2009h : subsampled stars outside this region are kinematic members of the stellar halo despite their spatial location 
near the surface brightness peak at the disc's mid-plane. In this hydrodynamical test case, red points mark original disc stars and purple 
points denote stars formed during the simulation's evolution. 



disc is less massive and thus more susceptible to the strong 
global resonance modes induced by the high-latitude accre- 
tion event. We note that in the hydrodynamical regime, the 
increase in friction accelerates the subhalo infall and thus 
results in less material stripped from the satellite galaxy at 
large heights above the primary galaxy plane; these accre- 
tion events are typically complete a few hundred Myr before 
their collisionless counterparts. 



3.1 Applicability to the Milky Way and M31 

Recent surveys indicate that the Galactic stellar halo may 
somewhat possess a dual nature similar to that observed 
in M31, being comprised of an outer spheroid of metal- 
poor stars having a slightly retrograde rotation with re- 
spect to the Milky Way disc, and an inner component with 
a flatter axial ratio and a mild prograde rotation, as well 
as stars more metal-r ich than the classic e xtended halo 
by a factor of three ilChiba fc Beers! |2000| ; ICarollo et all 
120071 ; iKinman et al]|2007l ). Although most of the Galactic 
halo is composed of an ancient stellar population, younger 
and richer stars have been observed there for decades 
dGreenstein fc SargentJl974l ; ISiegel et al.ll2009l ) ; a small frac- 
tion of these objects have peculiar velocities such that their 
calculated flight times from the disc are longer than their ap- 
parent lifetimes, and thus some debate exists over whether 
these sta rs formed in s itu in the low-density gas of the stel- 
lar halo ilKeenanlll992|), or are possibly undergoing binary 
rejuvenation i|Peretsll2009h . or were ejected from the disc 
via dyn amical interactions either internal to a multiple star 
system (|Blaauwl Il96ll ; iGualandris et all |2004| ) or invo lving 
the massive black hole at th e center of the Milky Way (|Hillsl 
1 19881 ; lYu fc Tremaindl2003T). O n a la rger interaction scale, 
the simulations of Abadi et al.l (|2009l ) have shown that dis- 
rupting dwarf galaxies may leave behind stars with velocities 



similar to or exceeding the escape speed of the primary sys- 
tem, which may cause some degree of angular collimation 
and travel-time commonality in a hypervelocity population, 
as observed in the fast-moving Galactic halo stars clustered 
in the Leo constellation. 

The majority of runaway disc stars, however, are not 
these hypervelocity objects but instead have moderately 
hot kinematics similar to those of the Milky Way's classi- 
cal halo, althoug h they tend t o retain significant rotational 
velocity as well ((Martin 2006 ), and can even be found far 
above the Galactic pl ane (jRamspeck et al.l l200ll . see also 
lAllen fc Kinmanll2004l and references therein). In a spectro- 
scopic survey of stars belonging to both th e thick disc and 
the stellar halo, iNissen fc Schuster! (|2009h report a num- 
ber of stars with halo-like speeds, and unexpectedly high 
a-enrichment, surmising that this could impl y either a dissi- 
patio nal collapse mechanism in the halo (as in lGratton et al.l 
120031 ) or the signature of an accretion event involving a large 
satellite galaxy. 

In the context of the Milky Way, we can distinguish 
dynamically between stars in the thin/thick disc struc- 
ture and stars that have been heated into orbits similar 
to those found in the stellar halo by employing the anal- 
ysis of Toomre energy space, in which the radial (U) and 
vertical (W) space velocities are summed in quadrature 
(T = [U 2 + W 2 ] 1 ^ 2 ) and diagramme d ag ainst the rotational 
velocity V (|Sandage fc Foutslll987l , see ICarollo et all 120071 
for a discussion of the solar neighborhood as viewed in this 
parameter space). In Figure [4] we extract disc stars from a 
thin annulus sweeping through the solar locality and cen- 
tered on the mid-plane, with a radial width of 100 pc and 
thickness of 300 pc, and plot the kinematic properties of 
this population. Despite being located spatially near the 
mid-plane of the post-accretion thin-thick disc structure, the 
number of stars with thick-disc kinematics has increased sig- 
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Figure 5. Upper: Stellar radial velocities (from the edge-on perspective) as a function of distance along the major axis, for subsamples 
of the collisionless (black circles) and hydrodynamical (red squares for old stars and purple triangles for newly-formed stars) events 
involving planar infall (0 = 0°). Center: Velocities have been corrected into the kinematic frame of the dis c, i.e., the mean local value 
of the rotation curve has been subtracted from each point for ease of comparison to the analysis of M31 bv lChapman et al. | j2006f ) (see 
£|3.1I) , in which the authors define boundaries in this phase space for each kinematic component. Note that the inner few kpc are excised 
in these panels to relieve overcrowding in the region where the mean radial velocity decreases sharply, and that these panels represent an 
appropriately-sized subsample of the particle distribution. Lower: The relative fraction of stars in the thick disc (solid) and stellar halo 
(dotted) components defined by their radial velocity ranges, as a function of distance along the major axis; the innermost disc is restored 
in this panel and the entire particle distribution is represented. 



nificantly when compared to the dynamical composition of 
the initial galaxy models, regardless of the satellite galaxy's 
orbital inclination during infall. Moreover, non-trivial num- 
bers of disc stars near the plane at the solar radius have 
been heated so dramatically by the accretion event that 
their three-dimensional velocities are as high as those ob- 
served in the Galactic stellar halo. By comparison, the kine- 
matics of the accreted satellite stars are consistent with a 
mixed population of thick-disc and stellar-halo stars, as de- 
termined in large part by the rotational speed of the ac- 
creted material, which lags the primary disc by an amount 
that correlates with increasing infall inclination angle (see 
iPurcell et al. 2009b, for an extended discussion of this be- 
havior as it pertains to the dynamical relationship between 
accreted stars and accreted dark matter). 

Though the results presented in Figure[4]are intriguing, 
it appears highly unlikely that the Milky Way has recently 
suffered such a dramatic accretion event as those we fidu- 
cially investigate here; not only is the Galactic disc far too 
cold and thin to hav e undergone much g lobal heating in 
the last several Gyr (|Purcell et al.|[2009al ), it is also quite 
deficient in stellar mass and angular momen tum when com- 
pared to a sample of local spiral galaxies (H ammer et al.l 



2007). These constraints indicate that the accretion history 
of the Milky Way has been substantially more quiescent 
than would be expected for a typical halo of virial mass 
M host ~ 10 12 M Q in the ACDM cosmology. On the other 
hand, it is plausible to suggest that minor events such as 
the ongoing accretion and disruption of the Sagittarius dwarf 
galaxy may also be capable of heating some small portion of 
disc stars to ejection, especially in the outer disc where the 
disc's self-gravity is less powerful; we may safely assume that 
the Galaxy has undergone some number of accretions in this 
regime, these events being truly ubiquitous in a cosmological 
context. 

At least one member of the Local Group may repre- 
sent more nearly the cosmological norm; the inner spheroid 
of M31 has been chemically identified as having a stel- 
lar population that is largely of inter mediate age and a s 
metal-rich as that galaxy's thick disc (jDurrell et al.l [20041 ). 
It may well be that M31, having had a more active re- 
cent accretion history than the Milky Wa y (as surmised by 
iFerguson et~ail l2002; Gil bert et ail 2009b. based on density 
and metallicity variations in halo substructure), has a stel- 
lar halo that has been significantly contaminated by stars 
originally belonging to the primary disc. Some indications 
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do exist that this is the case; uniform and young stellar pop- 
ulations have been shown to ex hibit strong rotation ou t to 
proj ected radius R DTO i > 40 kpc (jRichardson et al]l2008l . see 
also iFerguson et al.| [2005) . In addition, the Giant Southern 
Stream and associated stellar features are clearly the tidal 
remnants of an accretion event involving a fairly massive 
satellite galaxy; estimates of the progenitor's properties sug- 
gest that it was a large system traveling on a highly-eccentric 
orbit of nearly co-planar orientation with respect to the M31 
disc, and the surface brightness of the tidal stream strongly 
implies that the event occurred within the past few Gyr 
ijlbata et al.ll2004l : iFardal et alJlioOrj ; iFont et al.ll2006l ). Our 
result for the simulated infall with orbital inclination 8 — 0° 
predicts a substantial amount of disc-star ejection associ- 
ated with such an accretion event, an expectation that is 
consistent with the relatively bright inner halo of M31. 

Along these lines, in Figure [5] we present an analogous 
consideratio n to that found in the s tudy of metal-poor M31 
halo stars bv lChapman et al . (2006., here after C06), in which 
the authors utilize the kinematic model of llbata et all (|2005l ) 
in order to correct observed heliocentric radial velocities for 
the systemic and rotational motion of the M31 disc as well as 
the disc's inclination. Following C06, we present line-of-sight 
velocities across the major axis for our 8 — 0° runs from the 
edge-on perspective (upper panels) and subtract the rota- 
tional motion (middle panels). C06 used velocity cuts to 
kinematically represent the thin disc (|w r ad| < 100 km/s), 
thick disc (100 < |n ra d| < 160 km/s), and stellar halo 
(«rad > 160 km/s) of M31; we adopt the same cuts here for 
comparison, including only heated stars in the left panels 
and only accreted stars in the right panels. We obtain sig- 
nificant relative fractions in each component across the full 
length of the major axis, quite similar to the results of C06 
for M31. Across the full length of the extended discs major 
axis, approximately ~ 10% of heated disc stars have radial 
velocities consistent with the thick disc criterion set by C06 
for M31, and ~ 5% are consistent with the kinematic stellar 
halo definition. In the inner disc, where the rotation curve 
drops sharply and the bulge dominates the stellar density 
(although the bulge stars have been removed from this fig- 
ure), the dynamical contamination is even more pronounced: 
the thick disc and stellar halo represent ~ 15 — 25% each of 
the radial velocity population. Stellar material deposited by 
the infalling satellite also populates these regions; around 
~ 10 — 20% of these stars fall outside the thin-disc range in 
radial velocity. 



4 DISCUSSION 

Observations of diffuse light in the inner halos of external 
galaxies have tended to inspire efforts to reconcile the stel- 
lar mass of the spheroid with the properties of an inferred 
progenitor satellite galaxy. However, our experiments clearly 
demonstrate that an inner halo component may be primarily 
composed of stars ejected to large heights above the galaxy 
plane during the disc heating process. In our post-accretion 
remnant systems, we obtain surface brightness profiles typi- 
fied by a power-law index n = —2.5 along the minor axis as 
shown in Figure [21 only slightly flatter between 5 and 10 kpc 
than the pro j ected profile found to be a general feature by 
IZibetti et all l|2004h in their image-stacking exercise involv- 



ing over a thousand SDSS galaxies. In each of our simulation 
endstates, ejected material originally belonging to the pri- 
mary stellar disc accounts for ~ 40 — 80% of the luminosity 
found above 5 kpc from the plane, with the accreted sub- 
halo's stars being distributed over a much larger volume and 
thus contributing minimally to the surface brightness of the 
inner halo. 

This morphological contamination of the inner halo oc- 
curs simultaneously with the energizing process that also sig- 
nificantly contributes stellar material to the kinematically- 
defined thick disc and stellar halo. Just as in models of stel- 
lar halo formation, the scenario of thick disc assembly most 
likely involves multiple mechanisms, although there is ac- 
tive debate regarding the order of importance in these phe- 
nomena. The wide range of kinematic behavior displayed 
by low-mass thick disc galaxies may indicate, especially for 
cases in which signatures of counter-rotation appear, that 
these systems form primarily via the accretion of stripped 
stars during a minor mer ger, or from star formation a ssoci- 
ated with such a merger (jYoachim fc DalcantonlboOSl ). 

On the other hand, there is solid photometric evidence 
that more massive disc galaxies also contain thick compo- 
nents that are simply too faint to discern kinematically, es- 
pecially in the absence of significant rotational lag in any 
particular stellar population. This co-rotation is a particu- 
lar hallmark of low-latitude subhalo accretion, since in this 
case the satellite's stars will kinematically blend more effi- 
ciently into the host's rotation, and such uniform dynamics 
are also a general feature of thick discs formed largely by 
merger-related heating processes, which can achieve line-of- 
sight velocity dispersions similar to those found in surveys 
of these systems. In the case of the Milky Way, the thick 
component's faintness relative to the thin disc, as well as its 
measurably lagging rotation, likely indicate that it is not a 
product of heating that has occurred in the last several Gyr; 
observations have yet to determine whether the kinematic 
structure of the thin/thick Galactic disc is typical among 
systems with similar morphology and stellar content. 

Similarly, it remains to be shown self-consistently that 
a cold and thin Galaxy analogue can arise and persist to 
z = after the thick disc's formation at early times during 
the process of heating the proto-Galactic disc. However, we 
may still make a simple argument in favor of disc heating 
being largely responsible for the complex kinematic distri- 
butions in today's solar neighborhood. As alluded to above, 
the total mass of the disc under assault is an integral factor 
in any accretion event; according to first-order expectations, 
the strength of the stellar self-gravity will largely determine 
the extent to which any deposited orbital energy (or excited 
global resonance) heats the system. As an approximation, we 
would therefore expect more severe dynamical and morpho- 
logical damage to be inflicted by an accretion event involving 
a lighter disc, such as the Galactic progenitor at intermedi- 
ate redshift, leaving behind some fast-moving relics of this 
era in nearby stellar populations at the present day. A more 
detailed exploration regarding the role o f disc mass in the 
accret ion-heating process can be found in lKa zantzidi s~t al.l 
(2010, in preparation). We also note here that although pre- 
vious work has indicated that multiple minor mergers are 
unlikely to cause more significantly more structural and dy- 
namical dama ge than that inflicted b y the most massive 
satellite alone (|Kazantzidis et aTllio OSl , it may be that mix- 
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ing processes and resonance modes are enhanced by these 
multiple infalls; we defer pursuit and discussion of these is- 
sues for future work. 

It appears that the formation of inner galactic halos 
may be heavily enhanced by the destructive heating under- 
gone by a stellar disc during the infall of a large satellite 
galaxy. As we improve our understanding of stellar halo 
composition in nearby disc galaxies, we may begin to as- 
sess the degree to which these diffuse spheroids are al- 
loyed by intermediate-age stars that are unlikely to have 
been part of the primeval halo constructed long ago from 
the debris of faint and metal-poor dwarf satellites. Future 
Galactic surveys, such as SEGUE-2 and APOGEE of the 
third iteration in the Sloane Digital Sky Survey (SDSS-III; 
I Weinberg et~ai1l2007h . will map the outer disc and stellar 
halo of the Milky Way to an unprecedented degree of accu- 
racy, and M31 observations will continue to refine our un- 
derstanding of that galaxy's complex inner structure. These 
efforts will exquisitely define the kinematics, density struc- 
tures, and chemical compositions of the various populations 
in each component, and should easily be able to distinguish 
the relative properties of ancient, metal-poor halo members 
as compared to younger, more metal-rich stars that could 
have been heated to ejection from the disc during a massive 
accretion event. 

We would like to thank Larry Widrow and John Du- 
binski for kindly making available the software used to set 
up the initial galaxy models, and David Weinberg for useful 
discussions. We also thank Nicolas Martin for very help- 
ful suggestions during the referee process, which improved 
the manuscript greatly. CWP and JSB are supported by 
National Science Foundation (NSF) grants AST-0607377 
and AST-0507816, and the Center for Cosmology at UC 
Irvine. SK is supported by the Center for Cosmology and 
Astro-Particle Physics at The Ohio State University. The 
numerical simulations were performed on the IA-64 clus- 
ter at the San Diego Supercomputing Center, with ancil- 
lary experiments performed on the GreenPlanet cluster at 
UC Irvine. This work was also supported in part by an al- 
location of computing time from the Ohio Supercomputer 
Center | |http://www.osc.edu/| |. 



REFERENCES 

Abadi, M. G., Navarro, J. F., & Steinmetz, M. 2006, MN- 
RAS, 365, 747 

Abadi, M. G., Navarro, J. F., & Steinmetz, M. 2009, ApJL, 

691, L63 
Abe, F., et al. 1999, AJ, 118, 261 

Allen, C, & Kinman, T. 2004, Revista Mexicana de As- 
tronomia y Astrofisica Conference Series, 21, 121 

Bell, E. F., et al. 2008, ApJ, 680, 295 

Belokurov, V., et al. 2006, ApJ, 642, L137 

Benson, A. J. 2005, MNRAS, 358, 551 

Blaauw, A. 1961, Bull. Astron. Inst. Netherlands, 15, 265 

Boylan-Kolchin, M., Springel, V., White, S. D. M., & Jenk- 
ins, A. 2009, larXiv:0911.4484l 

Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2001, 
ApJ, 548, 33 

Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., 



Kravtsov, A. V., Klypin, A. A., Primack, J. R., & Dekel, 

A. 2001, MNRAS, 321, 559 
Bullock, J. S., & Johnston, K. V. 2005, ApJ, 635, 931 
Brook, C. B., Kawata, D., Gibson, B. K., & Flynn, C. 2004, 

MNRAS, 349, 52 
Brook, C. B., Kawata, D., Gibson, B. K., & Freeman, K. C. 

2004, ApJ, 612, 894 
Carollo, D., et al. 2007, Nature, 450, 1020 
Chapman, S. C, Ibata, R., Lewis, G. F., Ferguson, 

A. M. N., Irwin, M., McConnachie, A., & Tanvir, N. 2006, 

ApJ, 653, 255 
Chiba, M., & Beers, T. C. 2000, AJ, 119, 2843 
Clarkson, W., et al. 2008, ApJ, 684, 1110 
Conroy, C, & Wechsler, R. H. 2009, ApJ, 696, 620 
Cooper, A. P., et al. 2009, larXiv:0910.3211l 
Cote, P., Marzke, R. O., West, M. J., & Minniti, D. 2000, 

ApJ, 533, 869 
de Jong, R. S., et al. 2007, IAU Symposium, 241, 503 
de Jong, R. S., et al. 2007, ApJL, 667, L49 
De Lucia, G., & Helmi, A. 2008, MNRAS, 391, 14 
Durrell, P. R., Harris, W. E., & Pritchet, C. J. 2004, AJ, 

128, 260 

Fardal, M. A., Babul, A., Geehan, J. J., & Guhathakurta, 

P. 2006, MNRAS, 366, 1012 
Ferguson, A. M. N., Irwin, M. J., Ibata, R. A., Lewis, G. F., 

& Tanvir, N. R. 2002, AJ, 124, 1452 
Ferguson, A. M. N., Johnson, R. A., Faria, D. C, Irwin, 

M. J., Ibata, R. A., Johnston, K. V., Lewis, G. F., & 

Tanvir, N. R. 2005, ApJ, 622, L109 
Font, A. S., Johnston, K. V., Guhathakurta, P., Majewski, 

S. R., & Rich, R. M. 2006, AJ, 131, 1436 
Governato, F., Willman, B., Mayer, L., Brooks, A., Stin- 

son, G., Valenzuela, O., Wadsley, J., & Quinn, T. 2007, 

MNRAS, 374, 1479 
Gratton, R. G., Carretta, E., Desidera, S., Lucatello, S., 

Mazzei, P., & Barbieri, M. 2003, A&A, 406, 131 
Greenstein, J. L., & Sargent, A. I. 1974, ApJs, 28, 157 
Gualandris, A., Portegies Zwart, S., & Eggleton, P. P. 2004, 

MNRAS, 350, 615 
Gilbert, K. M., Font, A. S., Johnston, K. V., & 

Guhathakurta, P. 2009a, ApJ, 701, 776 
Gilbert, K. M., et al. 2009b, ApJ, 705, 1275 
Guhathakurta, P., Ostheimer, J. C, Gilbert, K. M., Rich, 

R. M., Majewski, S. R., Kalirai, J. S., Reitzel , D. B., & 

Patterson, R. J. 2005, |arXiv:astro-ph/0502366| 
Hammer, F., Puech, M., Chemin, L., Flores, H., & Lehnert, 

M. D. 2007, ApJ, 662, 322 
Helmi, A. 2008, A&A Rev., 15, 145 
Hills, J. G. 1988, Nature, 331, 687 

Ibata, R., Chapman, S., Ferguson, A. M. N., Irwin, M., 
Lewis, G., & McConnachie, A. 2004, MNRAS, 351, 117 

Ibata, R., Chapman, S., Ferguson, A. M. N., Lewis, G., 
Irwin, M., & Tanvir, N. 2005, ApJ, 634, 287 

Johnston, K. V., Hernquist, L., & Bolte, M. 1996, ApJ, 
465, 278 

Johnston, K. V., Bullock, J. S., Sharma, S., Font, A., 

Robertson, B. E., & Leitner, S. N. 2008, ApJ, 689, 936 
Kalirai, J. S., et al. 2006, ApJ, 648, 389 
Katz, N. 1992, ApJ, 391, 502 

Kazantzidis, S., Bullock, J. S., Zentner, A. R., Kravtsov, 

A. V., & Moustakas, L. A. 2008, ApJ, 688, 254 
Kazantzidis, S., Zentner, A. R., Kravtsov, A. V., Bullock, 



Heated Disc Stars in the Stellar Halo 9 



J. S., & Debattista, V. P. 2009, ApJ, 700, 1896 
Kazantzidis, S., Purcell, C. W., Bullock, J. S., Mayer, L., 

& Dotti, M. 2010, ApJ, in preparation 
Keenan, F. P. 1992, QJRAS, 33, 325 
Khochfar, S. & Burkert, A. 2006, A&A, 445, 403 
Kinman, T. D., Cacciari, C, Bragaglia, A., Buzzoni, A., & 

Spagna, A. 2007, MNRAS, 375, 1381 
Koch, A., et al. 2008, ApJ, 689, 958 

Kregel, M., van der Kruit, P. C, & de Grijs, R. 2002, MN- 
RAS, 334, 646 
Kuijken, K., & Dubinski, J. 1995, MNRAS, 277, 1341 
Launhardt, R., Zylka, R., & Mezger, P. G. 2002, A&A, 384, 
112 

Lequeux, J., Combes, F., Dantel-Fort, M., Cuillandre, J.- 

C, Fort, B., & Mellier, Y. 1998, A&A, 334, L9 
Maccio, A. V., Dutton, A. A., van den Bosch, F. C, Moore, 

B., Potter, D., & Stadel, J. 2007, MNRAS, 378, 55 
Martin, J. C. 2006, AJ, 131, 3047 

McConnachie, A. W., Chapman, S. C, Ibata, R. A., Fer- 
guson, A. M. N., Irwin, M. J., Lewis, G. F., Tanvir, N. R., 
& Martin, N. 2006, ApJ, 647, L25 

McConnachie, A. W., et al. 2009, Nature, 461, 66 

Minniti, D., Borissova, J., Rejkuba, M., Alves, D. R., Cook, 
K. H., & Freeman, K. C. 2003, Science, 301, 1508 

Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 
319 

Morrison, H. L., Miller, E. D., Harding, P., Stinebring, 

D. R., & Boroson, T. A. 1997, AJ, 113, 2061 

Moster, B. P., Maccio', A. V., Somerville, R. S., Johansson, 

P. H., & Naab, T. 2009. larXiv:0906.0764l 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 

462, 563 

Nissen, P. E., & Schuster, W. J. 2009, IAU Symposium, 
254, 103 

Perets, H. B. 2009, ApJ, 698, 1330 

Purcell, C. W., Bullock, J. S., & Zentner, A. R. 2007, ApJ, 
666, 20 

Purcell, C. W., Bullock, J. S., & Zentner, A. R. 2008, MN- 
RAS, 391, 550 

Purcell, C. W., Kazantzidis, S., & Bullock, J. S. 2009a, 

ApJL, 694, L98 
Purcell, C. W., Bullock, J. S., & Kaplinghat, M. 2009b, 

ApJ, 703, 2275 
Ramspeck, M., Heber, U., & Moehler, S. 2001, A&A, 378, 

907 

Read, J. L, Lake, G., Agertz, O., & Debattista, V. P. 2008, 

MNRAS, 389, 1041 
Rich, R. M., Reitzel, D. B., Howard, C. D., & Zhao, H. 

2007, ApJL, 658, L29 
Richardson, J. C, et al. 2008, AJ, 135, 1998 
Sackett, P. D., Morrison, H. L., Harding, P., & Boroson, 

T. A. 1994, Nature, 370, 441 
Sandage, A., & Fouts, G. 1987, AJ, 93, 74 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 
Siegel, M. H., Karata§, Y., & Reid, I. N. 2009, MNRAS, 

395, 1569 

Stadel, J. G. 2001, PhD thesis, University of Washington 

Starkenburg, E., et al. 2009, ApJ, 698, 567 

Stewart, K. R., Bullock, J. S., Wechsler, R. H., Mailer, 

A. H., & Zentner, A. R. 2008, ApJ, 683, 597 
Stinson, G., Seth, A., Katz, N., Wadsley, J., Governato, F., 

& Quinn, T. 2006, MNRAS, 373, 1074 



Stinson, G. S., Dalcanton, J. J., Quinn, T., Gogarten, S. M., 
Kaufmann, T., & Wadsley, J. 2009, MNRAS, 395, 1455 

van Zee, L., Barton, E. J., & Skillman, E. D. 2004, AJ, 128, 
2797 

Velazquez, H. & White, S. D. M. 1999, MNRAS, 304, 254 
Villalobos, A., & Helmi, A. 2008, MNRAS, 391, 1806 
Wadsley, J. W., Stadel, J., & Quinn, T. 2004, New Astron- 
omy, 9, 137 

Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, 
A. V., & Dekel, A. 2002, ApJ, 568, 52 

Weinberg, D. H., et al. 2007, Bulletin of the American As- 
tronomical Society, 38, 963 

Widrow, L. M., Pym, B., & Dubinski, J. 2008, ApJ, 679, 
1239 

Yoachim, P., & Dalcanton, J. J. 2006, AJ, 131, 226 
Yoachim, P., & Dalcanton, J. J. 2008, ApJ, 682, 1004 
Yu, Q., & Tremaine, S. 2003, ApJ, 599, 1129 
Zentner, A. R. 2007, International Journal of Modern 

Physics D, 16, 763 
Zibetti, S., & Ferguson, A. M. N. 2004, MNRAS, 352, L6 
Zibetti, S., White, S. D. M., & Brinkmann, J. 2004, MN- 
RAS, 347, 556 

Zolotov, A., Willman, B., Brooks, A. M., Governato, F., 
Brook, C. B., Hogg, D. W., Quinn, T., & Stinson, G. 
2009, ApJ, 702, 1058 



